Simultaneous cooling and entanglement of mechanical modes of a micromirror in an 

optical cavity 
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Laser cooling of a mechanical mode of a resonator by the radiation pressure of a detuned optical 
cavity mode has been recently demonstrated by various groups in different experimental configura- 
tions. Here we consider the effect of a second mechanical mode with a close, but different resonance 
frequency. We show that the nearby mechanical resonance is simultaneously cooled by the cavity 
field, provided that the difference between the two mechanical frequencies is not too small. When 
this frequency difference becomes smaller than the effective mechanical damping of the secondary 
mode, the two cooling processes interfere destructively and cavity cooling is suppressed in the limit 
of identical mechanical frequencies. We show that also the entanglement properties of the steady 
state of the tripartite system crucially depend upon on the difference between the two mechanical 
frequencies. If the latter is larger than the effective damping of the second mechanical mode, the 
state shows fully tripartite entanglement and each mechanical mode is entangled with the cavity 
mode. If instead the frequency difference is smaller, the steady state is a two-mode biseparable 
state, inseparable only when one splits the cavity mode from the two mechanical modes. In this 
latter case, the entanglement of each mechanical mode with the cavity mode is extremely fragile 
with respect to temperature. 

Keywords: Mechanical effects of light, Entanglement 



INTRODUCTION 



Mechanical resonators at the micro- and nano-meter scale are now widely employed in the high-sensitive detection 
of mass and forces (l|. Among the applications that have become possible are measurements of forces between 
individual biomolecules 0, forces arising from magnetic resonance of single spins and perturbations that arise 
C$ • from mass fluctuations involving single atoms and molecules [4|. The recent improvements in the nanofabrication 
techniques suggest that in the near future these devices will reach the regime in which their sensitivity will be limited 
by the ultimate quantum limits set by the Heisenberg principle, as first suggested in the context of the detection 
of gravitational waves by the pioneering work of Braginsky and coworkers [5[ . An important step in this direction 
would be the demonstration of cooling such microrcsonators to their quantum ground state. It would represent 
a remarkable signature of the quantum behavior of a macroscopic object, allowing to shed further light onto the 
quantum-classical boundary [f|. Recent experiments achieved promising results via cryogenic cooling via back- 
action cooling, in which the off-resonant operation of the cavity results in a retarded back action on the mechanical 
' system [§l 0, fiol [ill Il2[ [l3| , or by cold-damping quantum feedback where the oscillator position is measured through 
| a phase-sensitive detection of the cavity output and the resulting photocurrent is used for a real-time correction of the 
dynamics 14 , li|, 1(| 17 , 3| ■ As shown by recent theoretical results [lj| 2^, 21 , 2^, 2^, 24 , 2^, 2(| , these cooling scheme 
are in principle capable to bring the microresonator down to its ground state. Both mechanisms achieve cooling by 
increasing the damping of the mechanical resonator so that it becomes insensitive to thermal noise. It is important 
to analyze all the possible limitations to ground state cooling, also because this would open up the experimental 
realization of a number of genuinely quantum phenomena, such as the possibility of entangling an acoustic mode to a 
cavity quantum field [27j |. or to another mechanical oscillator 28, 

SHE HH, and even continuous variable quantum 
■ information protocols such as quantum teleportation [32| . and entanglement swapping [33j |. 

Here we shall consider only back-action cooling and the prototypal situation of an optical Fabry-Perot cavity with 
a rigid massive mirror at one end and a lighter, vibrating mirror at the opposite end. An intense laser beam drives 
a single, well separated, cavity mode which excites many internal vibrational modes of the oscillating mirror, via the 
radiation pressure. All the theoretical treatments of back-action cooling have up to now focused on a single cavity 
mode-single mechanical mode interaction, which is justified when a bandpass filter in the detection scheme is used, 
because coupling between the different vibrational modes is typically negligible. 
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In this paper we extend the treatment of cooling by including the effect of secondary acoustic modes whose resonance 
frequency is not far from that of the mechanical mode of interest, so that they cannot be neglected and filtered out by 
the detection process. We specialize to the case of a single additional mode and study its effect on cooling and on the 
entanglement properties of the steady state of the system. We find that cooling of the main mode crucially depends 
upon the difference between the two mechanical resonance frequencies. If this difference is larger than the effective 
damping of the secondary mode, cooling is not crucially affected by the presence of the adjacent mode, and the 
cavity mode is capable of simultaneously cooling the nearby mechanical mode close to its ground state. In this regime 
the cavity mode is entangled to each mechanical mode and the stationary state shows fully tripartite entanglement. 
If instead the two mechanical frequencies are closer than the effective mechanical damping, then the two cooling 
processes interfere destructively and each mechanical mode is no more cooled. This destructive interference affects 
also the stationary entanglement, which becomes extremely fragile with respect to temperature. The steady state of 



the system becomes a two- mode biseparable state, (i.e., of "class 3" 34|), which means inseparable only when one 
splits the cavity mode from the two mechanical modes. 

This paper is structured as follows. Section II introduces the physics of the optomechanical interaction inside an 
optical cavity and the linearized Langevin equations formalism. In Sec. Ill, we solve the dynamics of two mechanical 
modes coupled to the cavity field and compare it to the case of a single mechanical mode. In Sec. IV we characterize the 
simultaneous entanglement of the two acoustic modes with the cavity field and also study the tripartite entanglement 
of the steady state. Section V concludes the paper. 

SYSTEM DYNAMICS 

We consider an optical Fabry-Perot cavity of length L formed by a rigid massive mirror at one end and a vibrating 
micromechanical mirror at the opposite end, driven by a laser with frequency u)$. We shall refer from now on to this 
prototypal situation, even though the analysis could be easily adapted to other cavity geometries, such as the toroidal 
silica microcavities of Ref. 13; HI- The laser significantly drives only a single cavity mode with frequency uj c , from 
which it is detuned by Ao = lo c — llJq. The motion of the micro- mirror can be described by the set of its vibrational 
normal modes, each with its own resonance frequency Uj and damping rate jj. The Hamiltonian of the system is 

H = hw c a)a + J2 ^-(P 2 j + <fj) + H int + ihE^e-*""* - ae^ 4 ), (1) 
j 

where the cavity field annihilation operator a satisfies the commutation relation [a,a^] = 1, and the mechanical 
modes are described by dimensionless position and momentum operators satisfying [qk,Pj] = io~kj- Denoting by k the 
cavity decay rate, the parameter E is related to the input power P in by \E\ — y / 2Vi n n/hujo- The single cavity mode 
description is valid in the adiabatic limit when all the relevant mechanical frequencies uij are much smaller than the 
cavity free spectral range c/2L, which is typically satisfied for small cavities. In this limit the scattering of photons 
by the mirror motion from the driven mode to the other cavity mode is negligible [3a |. The interaction between the 
cavity mode and the vibrational modes is described by Hi nt and it is due to the radiation pressure acting on the 
surface S of the vibrating mirror. One has [3H ] 



S 



H int =- I d'rP^-uif), (2) 
where P(r) is the radiation pressure field and 



^1=Ey^-ft%(r) (3) 

is the displacement field of the mirror surface at point r. This field can be written as a sum over the corresponding 
(dimensionless) displacement field of each normal mode, Uj(r), which is characterized by an effective mass rrij = 
p J d 3 r \uj(r)\ 2 (p the mirror mass density). We consider a one-dimensional situation, i.e., we assume that the driving 
laser and the cavity are perfectly aligned. In this case, light is sensitive only to mirror surface deformations along the 
cavity axis, u x (r), so that Eq. ([2]) becomes 

H in t = ~ [ d 2 rP x (r)u x (r). (4) 



s 



In general, the radiation pressure due to an optical power V impinging on a mirror with reflection coefficient 1Z can 
be written as 



2V 

P x {r) = —Tlv 2 opt (r), 



(5) 



with v op t{r) denoting the spatial structure of the incident optical field on the mirror surface. Within the cavity, one 
can rewrite TP jc = h(u c /L)a^a and also assume 1Z ~ 1. One ends up with 



Hi n t = -h^2,G 3 Q a^aqj, 

3 

where the optomechanical couplings are given by 



(6) 



Qi = / n 



and 



(7) 



(8) 



is the overlap at the mirror surface between the cavity mode and the j-th mechanical mode. Due to the chosen 
normalization of v% pt (r) and Uj(r), the overlaps satisfy the condition < c,- < 1. Eqs. ©-© show that the radiation 
pressure directly couples the cavity mode only with the mirror collective displacement operator q e ff = G 3 qj. When 
the detection bandwidth involves only a single, isolated, vibrational normal mode of the microresonator, the collective 
coordinate q e f f is well approximated by the selected normal mode, and the single harmonic oscillator description 
usually adopted is justified. In the more general case, one has to include in the dynamical description of the system 
all the vibrational normal modes which contribute to the detected signal. 

The unavoidable action of damping and noise onto the dynamics associated with the Hamiltonian of Eq. |T]) is 
described by adopting the formalism of quantum Langevin equations 37i, 38] which, in the frame rotating at the laser 
frequency cl>q, ar e given by 



1j 

Pj 
a 



= -u^j ~ IjPj + G^ata + 
-(k + iA )a + iJ2j G J a aqj + E + \f2ka 



(9) 
(10) 
(11) 



The cavity input noise is delta correlated in the time domain (a m (t)a m '' (t')) — S(t — t'), while the mechanical 
Brownian stochastic forces with zero mean value t;j(t) are uncorrelated from each other and have the following, 
generally non-Markovian, correlation functions 



coth 



2k B T 



1 



(12) 



with ks the Boltzmann constant and T is the temperature of the reservoir of the micromechanical mirror. However, 
the involved mechanical frequencies are never larger than hundreds of MHz and therefore, as discussed in [25| (see 
also [39l]), even for cryogenic temperatures one can make the approximation 



7jt 



coth 



2k B T 



7r 



2k B T 



li ( 2u j + 1) • 



(13) 



where nj = [exp{huij / k B T} — 1] 1 is the mean thermal phonon number of mode j. As a consequence, the Brownian 
noise can be safely considered Markovian, that is, 



(2nj + l)S(t -t')+i 



S'(t - t') 



(14) 



where 5'(t — t') denotes the derivative of the Dirac delta. 



Ground state cooling is typically achieved when the radiation pressure coupling is strong. This can be obtained 
when the intracavity field is very intense, i.e., for high-finesse cavities and enough driving power. In this limit (and if 
the system is stable) the system is characterized by a semiclassical steady state with the cavity mode in a coherent 
state with amplitude a s (\ot s \ 3> 1), and a new equilibrium position for the vibrational modes, displaced by q 3 s . The 
parameters a s and q 3 s are the solutions of the nonlinear algebraic equations obtained by factorizing Eqs. (j9"|)- (|lip and 
setting the time derivatives to zero. They are given by 

qi = ^ (15) 

Pi = 0, (16) 
«- = ^Hfe> (17) 

where the effective detuning A(j\n is obtained from Ao by subtracting the frequency shift caused by the steady state 
radiation pressure 

A (A0 =A -|a s | 2 £^lI. (18) 

Ulj 

Then, we linearize Eqs. (j9l)- (fTT|) around the steady state values by writing operators as sums of averages plus 
fluctuations: a — a s + 5a, qj = q 3 s + 5qj and pj = p 3 s + Spj. The nonlinear terms Sa^Sa and SaSqj can be ignored when 
the fluctuations are much smaller than the mean value, and this is certainly satisfied when \a s \ S> 1. One therefore 
arrives at a system of linearized quantum Langevin equations 

5q =uj :j dp j , (19) 

Spj = —UjSqj — jjSpj + GjSX + £j, (20) 
SX = -kSX + A (N) SY + y/2KX in , (21) 

SY =~k6Y- A {N) SX + J2Gj6<lj + V2^Y m . (22) 

i 

We have chosen the phase reference of the cavity field so that a s is real and positive, we have defined the field 
quadratures SX = (5a + 8aJ)/>/2 and SY = (5a — 5aJ)/iv2 and the corresponding Hermitian input noise quadratures 
X m = (a m + a m '^j\[2 and Y m = (a m — a m '')/i^/2. We have also defined the effective optomechanical couplings 



^ mjLUjUo (k 2 + A 



r '"' V (23) 
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SIMULTANEOUS COOLING 



At the steady state, the energy of each mechanical mode can be written in terms of the variances of the corresponding 
position and momentum operators, 

U 3 = ^ [H) + (SP*)] • (24) 

In the absence of any cooling mechanism one has Uj = ftWj(rij + 1/2) and therefore when cooling is present we can 
write Uj = huij(n e ^ + 1/2), where n e ^ is the mean effective excitation number of the j mode, corresponding to an 
effective mode temperature 1^- = Fu+>j/[kB ln(l + )]• ^ one defines the vector of fluctuations 

u(t) = (5q 1 (t),5p 1 (t),...5q j (t),5pj(t)...,5X(t),5Y(t)) J , (25) 

and the vector of noises 

v(t) = (0,fi(i), ...0,€j(<)-. V2^6X m (t),V2^SY m (t)Y , (26) 
then Eqs. (fK))) - ({2"2")) can be written in a compact form as 

^-u(t)=A (N) u(t)+v(t), (27) 



where ^4(jv) is the drift matrix that governs the dynamics of the expectation values. Since the evolution is linear and 
the noise terms in Eqs. (fTT?|) - (|22[) are zero- mean quantum Gaussian noises, the steady state of the fluctuations is a 
zero-mean multipartite Gaussian state fully characterized by its correlation matrix V, whose elements are defined as 

(ui (oo)u m (oo) +u m (oo)u; (oo)) , . 

Vim = 2 • ( ' 

Using standard techniques [40j], one can determine the steady state correlation matrix V by solving the Lyapunov 
equation 

A {N) V + VA[ N) = -D, (29) 

where 

D = diag[0, 71 (2»i + 1) , ...0, 7j {2nj + 1) , «], (30) 

is the diagonal (2A + 2) x (2N + 2) diffusion matrix determined by the noise correlation functions. The stationary 
variances of the mechanical modes are given by the corresponding diagonal matrix elements of V. 
If only one mechanical mode is considered, the drift matrix assumes the following form 
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(31) 



where Am — Ao — G\/2u>\ is the effective single mode detuning. The stationary variances are given by (Sqf) = Vn 
and \Spi) — V22 and their exact expression is given in [25| . where they have been obtained by integrating the spectra 
obtained from the Fourier transform of the quantum Langevin equations. The two calculations coincide whenever 
a Markovian treatment of the quantum Brownian noise acting on the mechanical modes is made, i.e., when the 
approximation of Eq. Q13D is considered. 

As shown in 2^, 24 , |25|. |26| . cooling occurs when Am ~ u>i, i.e., when the laser is red-detuned with respect to the 
cavity mode, and the latter is resonant with the AntiStokes sideband of the laser. In fact, the laser light is scattered 
by the oscillating mirror into the Stokes and Antistokes sidebands with frequencies u>o ± u\. The generation of an 
AntiStokes photon takes away a vibrational phonon and is responsible for cooling, while the generation of a Stokes 
photon heats the mirror by producing an extra phonon. If the cavity is resonant with the Antistokes sideband, cooling 
prevails and one has a positive net laser cooling rate T given by the difference of the scattering rates, T = A- — A+. 
Refs. 23, 24, 25, 26| also show that laser cooling is optimized and can approach ground state cooling in the resolved 



band limit when k < ui± (actually when k ~ 0.2wi 23, 25|). Another important condition for ground state cooling is 
to have large optomechanical coupling G\, which is obtained for large intracavity power. However G\ has an upper 
bound imposed by the stability conditions [411. which in the case of a single mechanical mode and restricting to 
positive A(!), reduce to the single inequality 




7 (D = 1 _ 



G?A (a 



> 0. 



(32) 



However, as we have seen in Sec. II, the optical mode is always coupled to all the mechanical modes with a nonzero 
overlap Cj with the cavity mode at the mirror surface. Therefore, the actual stability conditions of the system are 
determined by the coupling with all the N excited mechanical modes, even the unobserved ones. By applying the 
Routh-Hurwith criterion [4l| . it is possible to see that, if we restrict to the cooling regime of positive detunings Any), 
there is always one nontrivial stability condition only, which is the direct A-mode generalization of Eq. 
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W = 1 _ 



A 2 



3 



G 2 

- S ->Q. 



(33) 



The violation of this condition leads to a bistable behavior, which has been experimentally verified in [4S 
assume that this stability condition is always satisfied from now on. 
In the case of two mechanical modes, the drift matrix is 
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(34) 



We have exactly solved the Lyapunov equation (|29p and analyzed the stationary position and momentum variances 
of the two mechanical modes in a parameter regime close to that of optimal cooling for a single mode, in order to 
see if and how the secondary mechanical mode affects ground state cooling. The results are shown in Fig. [U where 
the effective phonon number n e f f of the main mode (blue line) and of the secondary mode (green line) are plotted 
and compared to the single mode cooling case (red line). We have considered experimentally feasible parameters, (see 
caption) i.e., mechanical quality factors of the order of 10 5 and resonance frequency of the main mode uji/2it = 10 
MHz. 

We find two different situations, depending upon the value of the difference between the two mechanical frequencies, 
5u>2i = U32 — W\. Figs. [Ha) and[TJb) refer to the case when the two frequencies are well distinct, oj 2 — 1.7u>\ in (a) and 
U)2 = 2u>i in (b), and they plot n e ff versus the effective cavity detuning A( 2 ) = A — G\/2oj 1 — G 2 /2w 2 = Am — G 2 /2u; 2 . 
Fig.QJa) refers to the good cavity limit, k ~ 0.2wi, corresponding to a cavity finesse T — 1.5 x 10 5 , while Fig.^b) refers 
to a larger cavity bandwidth, k ~ u>i, corresponding to a finesse T = 3 x 10 4 . We consider a reservoir temperature 
T = 0.6 K, oscillators with mass m = 250 ng and a cavity length L = 1 mm. The results show that, when the two 
mechanical modes are well separated (S0J21 — u>i), the nearby mode does not disturb the cooling of the mechanical 
mode of interest, as witnessed by the perfect overlap of the blue and red curves, both in (a) and in (b). Even better, the 
secondary mode is simultaneously cooled close to its ground state (green curve). The comparison between Figs.[TJa) 
andQJb) shows that simultaneous cooling is influenced by the value of the cavity bandwidth k. In fact, the interval 



for the detuning A( 2 ) within which one has a significantly low value of n e f f is given by u>j — k < A 



can be seen from the width of the peak of the net laser cooling rate for mode j [23j, [24j, [25J, |2 
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k, as it 



(35) 



as a function of A( 2 )- Therefore, if the difference between the two mechanical frequencies is larger than k, the two 
modes are optimally cooled at two well distinct values of A( 2 ) and one can efficiently cool both modes only by fixing 
the detuning within a very narrow interval halfway between the two mechanical resonances, A/ 2 j ~ {u)\ + L02) /2 
(see Fig. (Ha)). However, due to the small value of k, the achievable value of n e ff at this intermediate detuning is 
appreciably larger (n e // — 0.5) than the optimal one achievable if one wanted to cool only one mode (n e // — 0.15). 
Instead, for a larger cavity bandwidth, k ~ ^w 2 i, a good simultaneous cooling of both modes is achievable in a 
significantly wider interval of detunings A( 2 ) (see Fig.[TJb)). 

Figs. QJc) and Hid) show that the situation is very different when the two mechanical resonances get very close. 
Fig. HJc) refer to w 2 = 0.95 x u>x, and one can see that in this case the occupancies n e ff of the two modes are 
appreciably higher than the one corresponding to a single isolated mode (red line). This means that when 5u>2i is 
small enough, the two cooling processes tend to interfere destructively. This is clearly confirmed by Fig. [TJd) , where 
n e ff at the fixed optimal detuning A( 2 ) = Wi is plotted versus the ratio oj 2 /u>i'- an occupancy of n e ff — 0.22 is 
reached for well separated frequencies but both modes are practically uncooled in a small interval around w 2 /cji ~ 1. 

This fact, at first sight unexpected, can be explained in terms of classical destructive interference between two 
resonant oscillators. A first explanation can be obtained by looking at the mechanical susceptibility of the main 
oscillator in the presence of the second mode, x* m (aj). It can be derived by Fourier-transforming the quantum 
Langevin equations, and it is given by 



[xT l M] _1 = [xi MP 1 - X2 M z 2 H g\gI 



where 



A 



A 2 ' 



(36) 



(37) 



\Xi (w)]" 1 = [Xo ( w )] 1 - ^ (w) G 2 , i = 1,2 



(38) 



(a) 



(b) 




FIG. 1: a) Simultaneous cooling of two acoustic modes (blue line: main mode, green line: secondary mode) of a single mirror 
versus normalized detuning A(2)/wi. The mechanical parameters are: 7i/27r = 72/27r = 100 Hz, lo\/2-k = 10 MHz, 102 = 1.7xwi , 
mi = 7712 = 250 ng and T = 0.6 K, that corresponds to initial occupancies m = 1250 and 712 = 735. The cavity of length L = 1 
mm and finesse T = 1.5 x 10 5 is driven by a laser of wavelength Ao = 1064 nm and power Pq = 30 mW (at A( 2 ) = u/i), which 
gives k ~ Gi — wi /5 . A red line is plotted but not visible owing to its almost complete overlap with the blue line, that portraits 
the behavior of the main mode cooling in the absence of the secondary mode. The final effective temperatures achieved are 
1 ~ 0.24 mK and — 0.42 mK. At the intersection of the cooling curves the effective temperatures are T*^ — 0.46 

mK and T^ !! — 0.79 mK b) Simultaneous cooling for a smaller cavity finesse T — 3 x 10 4 (corresponding to k ~ wi) and 
P = 100 mW (that gives Gi = 0.6 x wi). The main mode is chosen as before and its independent cooling curve is shown in 
red (again the two curves for 7i e // overlaps almost everywhere), while the secondary mode is fixed at frequency 102 = 2u>i- c) 
Simultaneous cooling for closely spaced modes u>2 = 0.95 x u>i and parameters n ~ wi/2, Gi ~ 0.3 x The cooling curve of 
the main mode (blue line) in the presence of the secondary mode (green line) is quite different from the independent cooling 
curve (red line), d) For the same parameters as in (c) the occupancy of the both modes at A( 2 ) = wi is shown versus 102 /uo\. 
An optimal occupancy of 0.22 is reached for well separated frequencies but it is strongly disturbed around 1^2/^1 — 1. When 
the two frequencies are equal cooling is practically absent. 



is the susceptibility of mode i modified by the radiation pressure of the cavity field, in the absence of the other mode, 
and 

[ ) C {u>)]- 1 = -[{^-^)-i^ i ] (39) 

is the bare susceptibility of the isolated microresonator. One can get an intuitive idea of the response of the mechan- 
ical mode of interest by rewriting (u) as the susceptibility of an harmonic oscillator with frequency-dependent 



resonance frequency (u) and damping 7 J e/J (w) 
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( w , e 'V) 2 -" 2 ) -iu^H 



(40) 



These two functions in the case of mode 1 arc plotted in Fig. [5] for the case of identical mechanical frequencies, 
u)% = u>2 ■ In this case both the effective damping and the effective frequency are strongly modified by the presence of 
the second mode (blue curve). The modification of the mechanical frequency due to radiation pressure is the so-called 



"optical spring effect" , which may lead to significant frequency shifts in the case of low-frequency oscillators [ill ] , 
but does not have significant effects in the case of higher frequencies, such as those of Refs. [g, [9|, llfj, and assumed 
here (see Fig. [UJb)). What is relevant is the modification of the effective damping that, in the presence of the second 
mode, quickly drops to a very small value in a narrow interval around resonance. Since cooling is signalled by an 
increased mechanical damping, this drop is just a manifestation of the suppression of cooling taking place when the 
two mechanical modes are resonant (see Fig. Et^)). This behavior is well described by the analytic expression of 
the effective frequency-dependent damping 7?^ (a;). The latter assumes a simple and transparent form when the 
susceptibility of mode i in the absence of the other mode can be taken as that of a usual resonator with an unmodified 
resonance frequency and a frequency- independent damping rate given by the net laser cooling rate T i: 



2 ) 



In fact, in this case, inserting Eq. ([4Tj) into Eq. (|36|) . one arrives at 
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(41) 



(42) 



Since the effective resonance frequency is not altered, the effective damping of mode 1 in the presence of the second 
mode is essentially determined by j^^(u!i), i.e., Eq. (|42|) evaluated at resonance. Therefore one has that when 
u>i = u>2, 7 4 e -^(wi) ~ 71 + 11(72/^) — 71 (since typically 72 <C T2), implying that at resonance the second order 
scattering processes mediating the interaction between the two mechanical modes completely suppress cooling. This 
process is a classical destructive interference phenomenon, similar to the classical analogue of electromagnetically 
induced transparency realized with two coupled oscillators in Ref. 43]. Another important information provided by 
Eq. ([4"2"|) is that this destructive interference starts to affect cooling when S1U21 < ^2, that is, Eq. (|4"2"|) shows that the 
"bandwidth" within which one has suppression of cooling is given by the effective mechanical damping, modified by 
the radiation pressure of the cavity, F2. 

One can see also view the suppression of cooling taking place when the two mechanical modes are resonant in a 
different way, starting from the interaction Hamiltonian of Eq. ©, which shows that the cavity mode directly interacts 
only with the collective coordinate Y2j ^olj- This suggests that, in the case of two mechanical modes, it is useful to 
pass to the effective "center-of-mass" and "relative" coordinates 
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With the new coordinates, the free Hamiltonian of the two mechanical modes becomes 
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where cu cm = {[GftV + [G 2 ] 2 ^} / {[Gl? + [G 2 ] 2 }, and u r = {[GI?lo 2 + [G 2 ] 2 Wl } / {[Gl? + [G 2 ] 2 }. This shows 
that at resonance, U2 = u>i, the relative coordinate q r is decoupled from the center-of-mass and therefore also from the 
cavity mode. As a consequence, it remains in its initial thermal state at the reservoir temperature and it is uncooled. 
Therefore, even though the center-of-mass is cooled close to its ground state, the two mechanical modes, 1 and 2, 
are only negligibly cooled because their steady state energy is determined by a weighted sum of the center-of-mass 
and relative mean energy. If instead u>x =/= LU2, the relative motion is coupled to the center-of-mass and the cooling of 
the latter by the cavity mode is able to partially cool also the relative motion. This "sympathetic" cooling is more 
efficient for larger coupling, i.e., for increasing 102 — ui, provided that both modes are not too far from resonance with 
the cavity. 
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FIG. 2: Plot of normalized effective mechanical damping rate (a) and mechanical frequency (b) of the main mode in the absence 
(red line) and presence (blue line) of a secondary mode as a function of u/uii. The parameters are yi/2n = 72/271 = 100 Hz, 
wi/27r = 10 MHz, ui2 = ui, mi = vti2 = 250 ng, T — 0.6 K, L = 0.5 mm, k = 0.3 x u}\ and Gi = wi/5. The quantities are 
plotted in the optimal cooling regime where Am = ui\ for the independent cooling case and A™ = u)\ for the simultaneous 
cooling case. At uj = u)x, the effect of the radiation pressure on the main mode is suppressed owing to the coupling to the 
secondary mode. 



ENTANGLEMENT PROPERTIES OF THE STEADY STATE OF THE SYSTEM 



From the results of Ref. [25} and [27| one can see that in the optimal cooling regime for a single mechanical mode 
A(x) ~ uji one has also a significant stationary entanglement between the mechanical and the optical cavity mode, 



which is also quite robust against temperature 27J. It is therefore interesting to study the entanglement properties 
of the stationary state of the tripartite system formed by the two close mechanical modes and the cavity mode. In 
particular, it is interesting to see if each mechanical mode is entangled with the optical mode, as in the single mode 
case, and also if the common interaction with the cavity mode enables to establish purely mechanical entanglement 
between the two vibrational modes. 

The 6x6 steady state correlation matrix V defined by Eq. I|28p. particularized for the case of two mechanical modes, 
can be written in terms of blocks of 2 x 2 matrices as 
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In order to quantify the bipartite entanglement of the Gaussian steady state of the three different instances of bipartite 
systems, we use the logarithmic negativity [3], defined as 

£ {Vup) = max{0, - In 2rT (V brp )}, (47) 

where Vu P is a generic 4x4 correlation matrix associated to the bipartite system of interest 

Vm p = (£ C B ) , (48) 



and 77 (Vbip) is given by [44 



v - (Vm p ) = (s (Vbip) ~ \/s (y Mp f - 4detV bip j ' , (49) 



with S (V^p) = det*4 + detS — 2detC The bipartite state is entangled if and only if (Vup) < 1/2, which is 
equivalent to the positive partial transpose (PPT) criterion, which is a necessary and sufficient entanglement criterion 
in the case of bipartite Gaussian states [4|| . 



To analyze the entanglement between one of the acoustic modes and the field, it suffices to eliminate the rows and 
columns that corresponds to the other mirror mode from the matrix V of Eq. (|46[) . We are left with a 4 x 4 matrix 



We have performed a numerical analysis around the parameter region considered in Ref. [2?| . which is within reach 
of state-of-the-art experiments (see the caption of Fig. [3]) and for which, in the presence of a single mechanical mode, 
one has a significant, stationary, optomechanical entanglement. 

We find that, similar to what happens for cooling, the entanglement properties of the steady state of the system 
strongly depend upon the value of the difference between the two mechanical resonance frequencies, 8u>2i- When 
the two modes are well separated, 8u)2i > Tj, the presence of the second mode does not affect too much the main 
mode-cavity field entanglement. Moreover, also the secondary mechanical mode is entangled with the cavity. This 
is illustrated in Fig. [3Ja) , where the optomechanical logarithmic negativity in the single mechanical mode case (red 
line) is plotted versus the cavity detuning and compared with the corresponding curves for the main mechanical mode 
(blue) and the secondary mode (green), when both modes are present. Optomechanical entanglement is more fragile 
than cooling because in the two mode case it is always smaller than that with a single mechanical mode. 

Figs. [3fb) and (c) show that the situation changes drastically when the two mechanical modes become very close 
in frequency, uj^/^x — 1. Both figures show the logarithmic negativity at a fixed detuning, versus the ratio W2/W1, 
at zero temperature (b) and at T = 0.4 K (c). We see that at zero temperature, the entanglement increases around 
resonance, 0J2/011 — 1. However, such an entanglement is very fragile with respect to temperature, and vanishes 
at T = 0.4 K (c) for a wide interval around the resonance condition. This behavior can be understood using the 
arguments of the preceding section. At mechanical resonance, the cavity mode is strongly coupled, and entangled, to 
the center-of-mass, and uncoupled from the relative coordinate. Modes 1 and 2 are linear combinations of q cm and 
q r and therefore their entanglement with the cavity mode is determined by both q cm and q r . At T = 0, both mode 
1 and mode 2 are entangled with the cavity mode thanks to the cavity-center-of-mass entanglement and because the 
quantum fluctuations of q r do not significantly affect it. However, as soon as temperature is increased, the thermal 
fluctuations of the uncoupled coordinate q r kill the entanglement of modes 1 and 2 with the cavity, even though 
cavity-center-of-mass entanglement is robust against temperature. Finally, the robustness of these optomechanical 
entanglements against temperature is analyzed in Fig. OJd). We see that even though never comparable to the case 
of a single mechanical mode, provided that the frequencies of the two modes are sufficiently far apart, i.e., 60121 <^ ^2 
(we have chosen L02 = 1.5o;i), one has a regime in which the two modes are simultaneously entangled, and this persists 
up to few Kelvins. 

One can also check if the two mechanical modes, even though not directly interacting, can become entangled at the 
steady state thanks to the common interaction with the cavity mode. Eliminating the entries in V that correspond 
to the cavity field, one is left with an all mechanical correlation matrix 

a)- <51) 

A numerical analysis of £ (V m _ m ) in a parameter regime around the region of optimal cooling, i.e., that of Figs. [T] 
and [21 shows no entanglement between the two mechanical modes, even when they are both strongly entangled to the 
same cavity field. Nonzero but extremely weak bipartite mechanical entanglement can be instead obtained in a regime 
where the oscillators are heavily damped and the cavity finesse is very high (« ~ 712 < wj. 2). This is consistent with 
the results of [46j |. where bipartite entanglement between two different macroscopic oscillators is analyzed. In fact, 
the present system is analogous to that of Ref. [3] , with the center-of-mass q r m of the two modes here playing the 
same role of the relative coordinate of Ref. f46j ]. As already shown in Ref. [46fl . purely mechanical entanglement is 
very fragile with respect to temperature and vanishes as soon as the occupancy of one of the modes is of the order of 
one. 



CLASSIFICATION OF TRIPARTITE ENTANGLEMENT 



The tripartite system under study can have various forms of tripartite entanglement and it is therefore important 
to classify the entanglement possessed by the steady state of the system. We determine the entanglement class of 
the system state by applying the results of Ref. 3J|, which has provided a necessary and sufficient criterion for 
the determination of the class in the case of tripartite CV Gaussian states and which is directly computable. This 




FIG. 3: Optomechanical entanglement, a) Mirror-field entanglement in the absence of secondary mechanical modes (red line) is 
plotted versus normalized detuning Am/ivi and compared with main mode-field entanglement (blue) and secondary mode-field 
entanglement (green). The parameters are 7i/27r = 72/27r = 100 Hz, U3i/2n — 10 MHz, W2 = 1.5 x u\, mi = m-z = 250 ng, 
T = 0.4 K, k = 0.9 x u)\ and Gi = u%. b) Enhancement of acousto-optical entanglement at zero temperature for close modes. 
The red line shows the value of the negativity of the main mode-field entanglement in the absence of the secondary mode, while 
the blue and green curves show the behaviour of the logarithmic negativity in the two-mode case when wi is fixed and u>2 is 
sweeped around ui. In this case Gi ~ 0.6o;i and we have fixed A2 = wi. c) The enhancement shown in b) is lost as soon as 
the temperature increases. Here the environment is at T — 0.4 K and the secondary mode frequency is varied between 0.5 x U\ 
and 3lji . d) Temperature robustness of entanglement in the collective case (blue and green lines for main and secondary mode 
respectively) compared to the independent case (red line). Here u>2 = 1.5 x uj\ and we have chosen A2 = L02, such that a 
simultaneous entanglement regime is obtained. 



classification criterion is mostly based on the nonpositive partial transposition (NPT) criterion proved in [47|, which 
is necessary and sufficient for 1 x JV bipartite CV Gaussian states. The NPT criterion of [13] can be expressed in 
terms of the symplectic matrix 

3 \ 

J = ©4 J < = (-lo) i=1 < 2 ' 3 ( 52 ) 
i=i ^ ' 

and of the partial transposition transformation acting on system k only. Transposition is equivalent to time 
reversal and therefore in phase space is equivalent to change the sign of the momentum operators. The NPT criterion 
states that a 1 x N CV Gaussian state is separable if and only if the "test matrix" Vk = AkVAk + ij /2 > 0. 
Therefore by evaluating the sign of the eigenvalues of the three possible matrices Vk, and using the NPT criterion, one 
can discriminate between the various entanglement classes. The results of this analysis are shown in Fig. [4] We find 
again that the properties of the steady state crucially depend upon the comparison between the difference between 
the two mechanical frequencies 5u2i and the effective mechanical damping Tj. When 5uj2i > Tj, the system is a fully 



tripartite state in a wide parameter region around the optimal cooling regime studied here, because the minimum 
eigenvalues for each bipartition are always negative (see Fig. UK a)). Instead, when the two frequencies are very close to 
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FIG. 4: Analysis of tripartite entanglement in the limit of well separated mechanical modes ui2 = 1.5 x W\ (a), and closely 
spaced modes U)2 = 1.01 x Wi (b). The minimum eigenvalues after partial transposition with respect to the main mode (red 
line), secondary mode (blue line) and field (green line) are plotted versus normalized detuning A( 2 )/o;i for a set of parameters 
7i/2tt = 72/2tt = 100 Hz, ^i/2tt = 10 MHz, mi = m 2 = 250 ng, T = 0.4 K, L = 0.5 mm, « = 0.5 x uii and Gi = 0.7 x u>i. 
For a large range of detunings, for well separated modes the total system is in a fully inseparable tripartite entangled state (a) 
while for close frequencies the resulting state is two-mode biseparable (b). 



CONCLUSIONS 



We have analyzed the effect of the radiation pressure of a cavity field mode on two vibrational modes of a mirror 
of Fabry-Perot cavity which are near in frequency and both within the detection bandwidth. We have considered 
the effect of the second mechanical mode both on ground state cooling via the back-action of the cavity mode, and 
on the optomechanical entanglement at the steady state. We have seen that the result crucially depends upon the 
difference between the two mechanical resonance frequencies. If this difference is larger than the effective bandwidth 
of the mechanical oscillators, given by the effective damping, the second mode not only does not affect cooling, but it 
is simultaneously cooled together with the main mode. Under the same condition, each mode is entangled with the 
cavity mode and the steady state is a fully tripartite entangled state. Instead when the two mechanical frequencies 
are very close to each other, cooling is destroyed by a classical interference effect and both modes are uncooled. In this 
condition, each mechanical mode is entangled with the cavity mode at zero temperature, but such an entanglement 
is extremely fragile with respect to temperature. Moreover, in the same regime, the steady state is a two-mode 
biseparable state which is inseparable only when the cavity mode is split from the rest. In these parameter regimes, 
the mechanical modes instead are never entangled; they can become entangled only for high finesse cavities, but the 
resulting entanglement vanishes at extremely low temperatures. 
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